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Recently the muonic hydrogen lamb shift has been measured with unprecedented accuracy, 
allowing for a precise determination of the proton radius. This determination is 5 sigma away 
from the previous COD ATA value obtained from (mainly) the hydrogen lamb shift and the 
electron-proton scattering. Within an effective field theory formalism, I will define the proton 
radius and briefly review some aspects of the theoretical prediction for the muonic hydrogen 
lamb shift, studying both the pure QED-like computation and the hadronic effects. 



1 Introduction 

The recent measurement of the muonic hydrogen Lamb shift, E(2P 5 / 2 (F = 2)) — E(2Sj/2(F : 
1)), 

E exp = 206.2949(32)meV 
and the associated determination of the electromagnetic proton radius [1]: 

(1) r p = 0.84184(67)fm . 

has led to a lot of controversy. The reason is that this number is 5 sigma away from the 
CODATA value, r p = 0.8768(69) fm [2]. If instead one uses this value in the theoretical 
expression of the muonic Hydrogen Lamb shift one gets the following discrepancy: 

(2) E exp - E th = 0.311 meV 

between theory and experiment. Two main options are clearly at hand: either the theoret- 
ical determination is not correct (or not as precise as claimed), or previous determinations 
of the proton radius were incorrect (or not as precise as claimed). Here we would like to 
study the theoretical expression of the muonic Hydrogen Lamb shift within an effective 
field theory perspective. We do it partially, and only focus on some few aspects, as a full 
analysis would require much more space. In particular spin effects will not be considered. 
We believe that the use of effective field theories helps in organizing the computation by 
providing with power counting rules to asses the importance of the different contributions. 
This will be even more important once higher order effects are included. For the present 
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discussion a 0(m r a 5 ) precision is enough to visualize the discrepancy. Higher order ef- 
fects are way smaller than the discrepancy found in Eq. (2). Moreover 0(m r a 5 ) is the only 
thing completely known at present 1 . 

The dynamics of the muonic hydrogen is characterized by several scales: 

m p ~ A x , 

m u m v 

n v ~ m n ~ m r — — — - 



m r CL ~ m e . 

By considering ratios between them the main expansion parameters are obtained: 

m n m p 1 

nip m v 9' 
m r cc m r tx. 2 1 



m r m r oi 137 

We use the effective field theory Potential Non-Relativistic QED (pNRQED) [3], Specially 
relevant for us is Ref. [4], which contains much more detailed information on the appli- 
cation of pNRQED to the muonic hydrogen, and we refer to it for details (see also [5-7]). 
pNRQED profits from the hierarchy m« 3> m }l 0i 3> m^a 1 and the Lagrangian reads 



P 2 

o 



2m r 



(3) L pNRQED = J d 3 rd 3 RdtS i (r, R, t) j id 

-V(r, V ,a lf r 2 )+er-K{R,t)^S{T,R,t)- J d 3 r^F fn ,F^ v , 

where S is the field representing the muonic hydrogen, R the center of mass coordinate 
and r the relative distance. 

V stands for the potential and admits an expansion in powers of 1 / m p : 

yM(r) V( 2 )(r) 



(4) F(r,p, £ r 1/£ r 2 ) = y(°)(r) 



m u ml 



The potential is obtained through matching to the underlying theory. Since pNRQED de- 
scribes degrees of freedom with E ~ wz ?i a 2 , any other degree of freedom with larger energy 
is integrated out. This implies treating the proton and muon in a non-relativistic fashion 
and integrating out pions. This is the step of going from Heavy Baryon Effective Theory 
(HBET) to NRQCD. By integrating out the scale m p a, pNRQED is obtained and the poten- 
tials appear. Schematically the path followed is the following: 

HBET(m n /triu) -» NRQED(ntua) -> pNRQED. 



1 This is also the precision presently reached in heavy quarkonium spectrum computations. This made 
provide with some cross checks between both systems. 



2 



XIV International Conference on Hadron Spectroscopy (hadron2011), 13-17 June 2011, Munich, Germany 



2 Pure QED contributions 

We first focus on the pure QED contributions. We mostly follow Pachucki's work [8-10], 
as it mainly follows an strict order by order in a computation, trying to accommodate their 
results in our formalism. See however [1] (or [11,12]) for more complete list of references. 

The static potential can be written in the following way in momentum space 

(5) y(°) = -±nZ }i Z p K V (k)^, 

(6) *'//W = * i + n(-k r 

where Tl(k 2 ) is the vacuum polarization due to electrons and can be computed order by 
order in a: 

n(k 2 ) = «n (1) (fc 2 ) + « 2 n (2) (fc 2 ) + a 3 n (3) (it 2 ) + ... 

(7) cc v (k) = cc eff (k) + £ Z^^fik) = cc eff (k) + Scc(k) , Sx(k) = O(^). 

n,m=Q 
n+m=even>0 

Ex, p E[, p' 

» i » 

i 
i 
i 

k , k 

i 
i 



Figure 1 : Leading correction to the Coulomb potential due to the electron vacuum polarization. 
k = p — p' and fcg = E\ — E\. 

The leading order contribution to the lamb shift comes from the one-loop vacuum polar- 
ization correction to the static potential (see Fig. 1) 

E LO = (n\SV\n) = 205.0074 mev = 0(m r a 3 ) . 

The C(m r a 4 ) contribution to the lamb shift comes from the two-loop static potential and 
from the iteration of the one-loop potential in quantum mechanics perturbation theory. 
The latter yields AE = 0.151 meV. The former is purely due to vacuum polarization 
corrections (see Fig. 2) and yields AE = 1.5079 meV. 
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Figure 2: The static potential at two loops. 

The three-loop static potential contribution due to the vacuum polarization (and the asso- 
ciated iterations from perturbation theory were computed in Ref. [13] (see also [14] for an 
small correction). The result was quite small AE = 0.0076 meV = 0(m r ct 5 ). 

All previous contributions were due to the vacuum polarization. The only contribution 
from the static potential that is not due to the vacuum polarization at O(m r 0i 5 ) comes from 
5a.. It is a light-by-light (Wichmann-Kroll and Delbriick) contribution and very small [15] 



AE ~ -0.0009 meV. 



It should be mentioned that the limit m e — > of the static potential is known at three loops 
from QCD [16, 17], which could be used as a check. It is also reasonable to think that the 
result with finite m e could also be obtained from these results (albeit numerically) with a 
finite amount of work. 

There are no corrections due to the l/m^ potential at O(m r oi 5 ). From the potential 
(see [4] for its expression in pNRQED) there are the tree level relativistic corrections, which 
give AE = 0.0575 meV = 0(ma 4 ). The incorporation of the one-loop vacuum polarization 
to the relativistic l/m\ tree-level potential gives the following result AE = 0.0169 meV = 
0{ma 5 ) [10]. 

In order to complete the pure QED O(moi 5 ) corrections one has to include the interaction 
with the ultrasoft photons (see Fig. 3). They yield the result (taken from [8]) 

AE = -0.6677 meV = 0{ma 5 ) . 

The -J- ultrasoft effects contribute 

Trl 

AE = -0.045 meV = 0(ma 5 ^-) . 

m v ' 

In pNRQED these results would not come from the interaction with the ultrasoft photons 
only, as it would be factorization scale dependent, they also include effects due to the 
NRQED matching coefficients encoded in the 1/tni potentials. The procedure is pretty 
much the same the one used for positronium in Ref. [5]. The details for muonic hydrogen 
will be worked out elsewhere. 
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Figure 3: self-energy correction to the muonic hydrogen energy due to the interaction with ultra- 
soft photons. 



At O(moc 5 —^) one starts to have overlap with hadronic effects, which we discuss in the 

Ttlp 

next section. 



3 Hadronic Contributions 

In the previous section we have considered the proton to be point-like. We now incorporate 
the finite-size effects due to the hadronic structure of the proton. These effects are encoded 
in the coefficient multiplying the delta potential (note that the combination of NRQED 
matching coefficients that appears in the potential is always the same). 

(8) < 2 j(r) = Id^W - AE = 1 DWI(^)3, M 



where 



(9) D h d ad = -c\ ad - 16nad h 2 ad + ^c^ rf . 

We define c 3 = c f nt ~ me + cf d , d 2 = d f nt - like + df d c D = ^'"^ + c ^d. so that c had > 
d\ ad , c^ d are the left-over of the matching coefficients of NRQED Lagrangian 

(10) 5C = ... %hv^ v + ■■■- e^N^V ■ ENp + • ■ ■ + 

after subtraction of the point-like contributions. We do in this way because traditionally 
the point-like contributions are already included in the "pure" QED corrections described 
in the previous section 2 . I more extended discussion can be found in Refs. [4,6]. 

flhad enco d es the hadronic vacuum polarization effect. Its contribution to the Lamb shift is 
tiny, AE = 0.011 meV, and not much subject to uncertainty as it can be determined with 
enough precision from dispersion relations. 



2 Note though, that for an strict effective theory point of view, at scales of the order of trip, it is not a good 
approximation to consider the proton point like. Therefore, in a way, we are introducing an "spurious" contri- 
bution in the hadronic matching coefficients. 
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More subject to discussion are the hadronic corrections associated to Cg . They are usu- 
ally split into two terms (see the discussion in Refs. [4,6]): c\ ad = c\ a ^ emach + c ! ^ o} . We 
symbolically draw them in Figs. 4 and 5, and discuss them in the next subsections. A 
common feature of both of them is that they are power-like chiral enhanced: ~ -f-. This 
is very important, as it allows chiral perturbation theory to predict the leading order term 
without introducing any extra parameter. The resulting correction to the Lamb shift is of 




Figure 4: Symbolic representation (plus permutations) of the Zemach (r 5 ) correction, Eq. (11). 




Figure 5: Symbolic representation ofEq. (17). 



3.1 Zemach correction, (r 3 ) 

It is the one analogous to the Zemach correction defined in the hyperfine splitting [18]. It 
is also common to rewrite it in terms of a coefficient (fp) 

rm r had --a. 2 m 2 m(r 3 ) ( r P _ 96 r D _ x 1 c (o) c (2) 

y ll > c 3,Zemach ~ ^ a m p m }'\ r p) > ^3 ~ ~ J a fc6 U E °E ' 
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where is the electric Sachs form factor to order n in the chiral counting (G^ = 1). 
We also use dimensional regularization (D = 4 + 2e). This gets rid of power-like diver- 
gences which are then automatically set to zero (no need for counterterms as one would 
with cutoff regularization). The final result is finite and it is possible to obtain an analytic 
expression for the leading term in the chiral and large N c expansion (by including the A 
particle contribution). It reads [4] 

+ i^T ; g (x) + &™ E (x 

where (A = M A - M p ~ 300 MeV) 

mi _ (-l)T(-3/2) f 2(r + 2) 1 

(13) Q " r(r + i)r(-3/2-r) \ B(,+2r ~ TUT* 4 **! ' r ~ °' 



f 2-n 







_ w!(2w-l)!!r[-3/2] 
" ~ 2(2n)!!r[l/2 + n] 



This expression produces the following number for (r 3 ,) and the associated energy shift: 

( r n)\xPT (T%) 

(15) v \ = 1.9 (Pineda) -> AE = 0.010^- = 0.019 meV 

fm 3 fm J 

This number can be compared with some recent determinations of (r 3 ) using dispersion 
relations [19-22] 



{ryiVexp" _ f 2.71(13) Friar -Sick 



fm 3 \ 2.85(8) Bernauer — Arrington 



AE = 0.027 - 0.029 



In principle the difference between these two determinations comes from different fit func- 
tions and data, which may give a first estimate of the associated uncertainty of the dis- 
persion relation analysis. We find quite reassuring that the difference with the chiral com- 
putation is around 40 %, which could be easily accommodated with higher order correc- 
tions. Much more difficult to accommodate would be the value advocated in Ref. [23], 
(r 3 ) ~ 36.5, from a direct fit to the muonic hydrogen Lamb shift using the COD ATA value 
for the proton value. This would require that higher order corrections in the chiral compu- 
tation to be a factor 15 larger than the leading order result. We believe this is at odds with 
chiral symmetry, even more so taking into account that one of the motivations of such pro- 
posal was the lack of experimental data at low momentum, but it is precisely in this region 
where chiral perturbation theory should work better. 
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3.2 Polarizability correction 



The determination of the polarizability correction from experiment is on more shaky grounds 
than for the Zemach correction, producing the larger uncertainty in the theoretical expres- 
sion for the Lamb shift. The reason is that dispersion relations do not fix the result com- 
pletely. The final number used in [1] was taken from the average in Ref. [24] AE = 0.015 ± 
0.004 using the results from, [9] AE = 0.012 ± 0.002 meV, [25] AE = 0.017 ± 0.004 meV, 
and [26] AE = 0.016 meV. For a recent discussion see Ref. [27]. 

Here again chiral computations may turn out to be crucial to asses the size of this cor- 
rection. The reason, as before, is that the polarizability correction is power-like chiral en- 
hanced. Therefore, chiral perturbation theory can predict the leading order term with no 
new parameter. This is the attitude followed in Ref. [7], where a chiral computation using 
dispersion relations yielded 



had 
L 3,pol 

where 



-e I m p m }l 



V" = iJi'xe"l'(p,s\TJ"(x)r(0)\p,s), 
which has the following structure (p = q ■ p/m): 



qVq v 



Si(p, q 2 



m- 



m 



2 1 



til 



VP 



q v )S 2 (p,q 2 



(16) 



in, 



e^q^A^p.q 2 



in- 



?J lv P a q p ((m p p)s a - (q • s)p a ) A 2 (p, q 2 ) 



After introducing the chiral expressions for the structure factors from the diagrams in Fig. 
6, one obtains 



(17) 



Jiad 



4 2 m V ( gA 

-e mi — - 1 

F m 7 



dw 



(In) - 1 (1 + k 2 ) 4 J 71 



w 



D-5 



x {(2+ (l + k 2 ) 2 )A E (a; 2 ,k 2 ) + (1 + k 2 ) 2 k 2 w 2 B E (a; 2 ,k 2 )} 
where (for D = 4) 



2 i 4 m /< l 

(l+k2)2 



(18) 



1 

47T 



1 

+ \/l + w 2 + / dx — = ~ X 

J + X 2 W 2 + X(l 



:(1 — x)w 2 k. 2 



8n 



dx 



2x 



y/1 + X 2 W 2 



x(l 



X)W 



2 k 2 



1 

- / dx — 
2 1 (1 



{l-x)(l-2x) : 



+ x 2 w 2 + x(l — x)w 2 'k 2 )^ 
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(Seagull) 



Figure 6: Diagrams contributing to T'L Crossed diagrams are not explicitly shown but 
calculated. 



This gives the number 
(19) 



AE 



X PT 



(pions) 



0.018 meV . 



We consider a more through chiral study of this object, in particular including the A parti- 
cle, compulsory. The introduction of the A particle produced a large effect in the case of the 
Zemach correction, something similar may happen here. Whereas we can (and should, see 
Ref. [28] for a recent discussion) further analyze the error associated to the polarizability 
correction, we would like to emphasize that in order to explain Eq. (2), the corrections to 
the leading order chiral computation should be a factor 15 larger than the number obtained 
in Eq. (19). 



3.3 Definition of the proton radius 

From the effective theory point of view, the proton radius corresponds to an specific combi- 
nation of the Wilson coefficients of the effective theory. Let us see how this relation appears. 
One first considers the following matrix element 



(20) 



(p',s\Ji l \p,s)=u(p') 



u{p). 



We are interested in the low energy limit of the form factors 

„2 



(21) 



W) = F i 



111 



2 t ' "• 



9 



XIV International Conference on Hadron Spectroscopy (hadron2011), 13-17 June 2011, Munich, Germany 



or more precisely of the Sachs form factors 

(22) G E (q 2 ) = F l {q 1 ) + ^U?ft, G M (q 2 ) = F 1 {q 1 ) + F 2 (q 2 ) . 

The proton radius is usually defined as the derivative of the Electric Sachs form factor at 
zero momentum: 



dG P/E (q 2 



(24) c D = 1 + 2F 2 + 8F{ = l + 8m 2 



dq 2 



q 2 =0 



This set of equations allows to visualize the relation between the proton radius and the 
matching coefficients of the effective theory [4]. They also highlight the problem of defin- 
ing the proton radius through the derivative of the Electric Sachs form factor at zero mo- 
mentum, as this object is infrared divergent. Then, what is the definition of the proton 
radius used to obtain Eq. (1)? It corresponds to subtract the point-like effect of c^: 



3 1 

r 2 v = j— (c D {v) - c DfP0int -uke{y)) 

4 TTl-p 



where, in the MS scheme, 



C D , point-like — 1 + — I ^ ^ n 



4 Conclusions 

We have briefly reviewed the theoretical determination of the spin-independent correc- 
tions to the Lamb shift that contribute at O(m r oc 5 ). We believe that it is important to have 
a model independent and efficient approach to the problem. Effective Field Theories are 
suitable for this task. The use of effective theories highlights that the proton radius is a 
matching coefficient of the effective theory and, in general, an scheme / scale dependent 
object. In principle, this is not a problem as far as one knows the definition one is using. 

The pure QED computation appears to be solid, not to say extremely reliable. Yet it could 
be interesting to reanalyze some parts of the computation from an effective theory perspec- 
tive. We also remark that the correction coming from the three-loop static potential has 
only been computed by one group. On the other hand, the analogous QCD static potential 
has been obtained by two groups. We believe that with a reasonable amount of effort, such 
computations may yield cross checks of the QED computation for muonic hydrogen. 
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For the hadronic corrections there are precise determinations of the effects associated to the 
Zemach correction, (r 3 ), using dispersion relations. Chiral perturbation theory provides 
with a highly non-trivial double check of the magnitude of this correction. The reason is 
that the chiral computation is power-like chiral enhanced. It actually linearly diverges in 
the chiral limit. Therefore, the leading order computation in chiral perturbation theory 
is a pure prediction, with no free parameter. This rules out much larger values of (r 3 ) 
than those obtained from experiment, as such values would be in tension with the chiral 
perturbation theory prediction. 

The polarizability correction is the major source of uncertainty. The reason is that dis- 
persion relations alone are not able to fully determine this quantity, suffering from some 
ambiguity in the parameterization. Therefore, the chiral perturbation theory result may 
turn out to be crucial here to determine the size of this correction. Again the chiral com- 
putation is power-like chiral enhanced and linearly diverges in the chiral limit. Thus, the 
leading order computation in chiral perturbation theory is a pure prediction, with no free 
parameter. At present there is room for improvement over the result obtained in Ref. [7] 
using chiral perturbation theory with dispersion relations. In particular, it does not include 
the contribution due to the A particle, which, in the case of the Zemach term, turned out to 
be quite important. It will then be very important to compute it to really asses the size of 
this correction. Obviously any eventual determination from lattice of this quantity would 
be most welcome. 
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